Inhomogeneous states in non-equilibrium superconductor/normal metal tunnel 
structures: a LOFF-like phase for non-magnetic systems. 
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We analyze non-equilibrium states in a tunnel superconductor- normal metal (N/S/N) structure 
in the presence of a tunnel current /. We use an approximation of an effective temperature T and 
CNj , calculate the current-voltage I-V characteristics. It is shown that the I-V dependence may have an S- 

shaped form. We determine nonuniform current I(x) and temperature T{x) distributions that arise 
as a result of instability of the uniform state with negative differential conductance (dl/dV < 0). We 
discuss an analogy with equilibrium superconductors with an exchange field in which nonuniform 
■ states predicted by Larkin/Ovchinnikov and Fulde/Ferrell are possible. 

; PACS numbers: 74.40.+k, 74.78.Fk, 74.25.Fy,72.20.Ht 
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^ ■ I. INTRODUCTION 

Study of superconductor/normal metal (S/N) or superconductor/superconductor hetero-structures has been a pop- 
ular topic during last decades pi 0, 0, 0] • Both the equilibrium and non-equilibrium properties of such systems were 
under intensive investigation. Non-equilibrium phenomena are of special interest because they are essentially differ- 
ent from the equilibrium ones but several interesting effects have been observed even in the cases when a deviation 
from the equilibrium is small. As an example, one can mention a non-monotonic behavior of the resistance of S/N 
structures as a function of temperature T or applied voltage V @, IS 0j fl 0] • 

If the applied voltage is not too small, the superconductor S (for example, in a N/S/N system) goes out of the 
' equilibrium and new phenomena come into play. Study of non-equilibrium effects in superconducting hetero-structures 
as well as in homogeneous superconductors in the presence of optical or microwaves radiations has been carried out 
during a long time (see reviews [l(| HH, E3] and references therein). For example, Eliashberg suggested to stimulate 
superconductivity by a microwave radiation that leads to an essential deviation of the distribution function / from 
the equilibrium one [3. As the energy gap A in S is related to the function / via the self-consistancy equation, the 
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ly^ I energy gap and the critical temperature may increase in the presence of a microwave radiation, which results in an 



X 



increase of the critical current [14j . Another interesting effect caused by a non-equilibrium distribution function is the 
sign reversal of the Josephson critical current j c in a multi-terminal SNS junction [151 ]. It turns out that, provided 
. the distribution function in the normal metal N controlled by additional N' electrodes differs significantly from the 
1 equilibrium one, the critical current j c becomes negative (estate) [l6l [l7l [l8j. 
0^ . In recent years, the interest in studies of non-equilibrium effects in S/N or S/S' hetero-structures has revived. This 
■ is due to the progress in nanotechnology, the possible applications of such structures in low-temperature devices fl9j 
and to the progress in theoretical research p], 0, 0, 0, H(| ■ 

Perhaps, the simplest system in which one can study non-equilibrium effects is a N/S/N structure with a bias 
current 7^. Such a structure was studied theoretically in recent papers [U H2]. In Ref.[2l| it was assumed that 
the resistance of the S/N interfaces is negligible. The calculated I-V characteristics (CVC) was shown to be of the 
so-called N-shape type, that is, three values of voltages correspond to one value of the bia s current i&. Earlier this 
type of the I-V characteristics was studied in other superconducting systems 
The opposite case of a large S /N interface resistance was considered in Ref . 



12 , 123|, 12J, [2| 

22J . In this case the I-V curve is of the 



S-type, i.e. three values of the current correspond to one value of the voltage V between the N leads. Such kind of the 
CVC was known from previous studies of the superconducting structures. For example, it can be realized in S/I/S 
junctions [T3, EH HH, H3 and in granular superconductors [28, 29]. The authors of Ref. [22j| investigated the stability of 
this system and came to the conclusion that the system is stable. However, this conclusion is valid only in the case of 
small lateral dimensions of the structure. From a general theory of the system with negative differential conductance 
(Gd = dl/dV < 0), it is known that the states with negative Gd are unstable and in case of the S-shape I-V curve 
a stratification of the current density occurs as a result of the instability [3(J E| • Inhomogeneous states in different 
superconducting systems were studied in approximate models with account for electron-electron or electron-phonon 
inelastic scattering [l(| [ll| . 
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In the present paper, we analyze the N/S/N structure of the type considered in Ref.[22j assuming however that 
lateral dimensions are large. The latter assumption turns out to be very important. We find in this limit that a new 
inhomogeneous state is possible, which contrasts the situation in small sized systems. Physics of this inhomogeneous 
state is very similar to that of the inhomogeneous state in equilibrium superconductors with an exchange field (the 
Larkin-Ovchinnikov-Fulde-Ferell states [33l. l34j). It is worth emphasizing though that, in the system considered here, 
there are no magnetic interactions and their role is played by a finite voltage V (see below). Under certain assumptions 
we calculate the I-V curve and study possible inhomogeneous states. Most of the results are obtained in an analytical 
form. In particular, we find the spatial distribution of the current density j(x) and energy gap A(x) and draw an 
analogy with those in the LOFF states. 



II. MODEL. BASIC EQUATIONS 

In order to calculate physical quantities we use a kinetic equation for the non-equilibrium distribution functions 
f±. The procedure of the derivation is quite standard (see, e.g. (36|) and is based on using quasi-classical Keldysh 
matrix Green functions g. They are expressed in terms of the two distribution functions f± as g = g (fa/- +fo/+) — 
(73/- + to/+).9" 4 , where T3 is the f z Pauli matrix and fo is the unit matrix. The functions /+ and /_ correspond to 
the functions Jl and fx introduced by Schmid and Schon [32J. 

The considered N/S/N system is shown in Fig. la. We assume that the tunnelling probability through the interface 
is small and that the retarded (advanced) Green's functions g R ( A > have the same form as in a bulk superconductor: 
g-R(A) _ g R(A)^ + fR(A) if2 with j-R(A) = A/ v /( e±i7 )2_ A 2 and g R(A) = - (pWf. The only difference is 
that, in contrast to the bulk case, the damping 7 caused by a weak proximity effect, enters the formula for the 
function f R ( A \ In this limit the damping 7 is equal to 7 ~ eo, where eo = D/(RsN&d), D = v 2 t/3 is the diffusion 
coefficient, Rsn is the SN interface resistance per unit area, a and d are the conductivity (in the normal state) and 
the thickness of the superconductor. The thickness of the S layer is assumed to be small (d « £s ~ y/D/A), so that 
all quantities do not depend on the z— coordinate. The opposite case was considered in Ref. [2l|, |4l|. The thickness 
of the N leads is supposed to be large (d/v >> £,s) so that the N metals are in the equilibrium state. 

The function /_ (e, V) appears due to a branch imbalance and, being dependent on the electric potential V, 
determines the current /. In case of a fully symmetric N/S/N system /_ = 0. The function /+ (e, V) determines the 
energy gap A. Since the function / + (e, V) depends on the applied voltage V, the gap A is also a function of V. The 
dependence A(V) is given by the self-consistency equation 



1 = XRe f D deAMl (1) 
Jo Ve 2 - A 2 ^ ' 

where e = (e + ij), X and 9d are the coupling constant and Debye energy, respectively. Eq. {T]) is a generalization of 
the conventional BCS equation to the non-equilibrium case. 

The functions f± can be found from the kinetic equation that can be written in S in the form (see for example [42] ) 



d(Ms(e)f±) d(D ± (e)df±) 



dt 



dx 2 



= ^ e aA a ±(e) + S in {f±} 



(2) 



a—l,r 



where 



•A/s(e) = Re{g R (e)} = Re 



v 7 ? 2 - A 2 



(3) 



is the normalized density-of-states (DOS), ti tT — (D / Ri yr ad) with Ri >r being the resistance of the left (right) S/N 
interfaces, and T>±(e) is an energy-dependent diffusion coefficient, 



D 



ee* T A 2 



v 7 ? 2 - AVe* 2 - A 2 



(4) 



with e* = e — ij. 

The functions Ai, r ±(t) in Eq. Q are defined as: 



A l!r ±(e)=Afs(e)[F ± (e,V lir )-f ± (e)] 



(5) 



3 



a) 





FIG. 1: a) The system under consideration, b) The proposed system in which an inhomogeneous distribution of the energy 
gap can be probed with a point contact spectroscopy. 



The distribution functions in the normal metals are assumed to have the equilibrium forms shifted by the applied 
voltages V. In a fully symmetric system, the voltages in the right (V r ) and left (V/) N metals are equal to V r — —Vi = V. 
In this situation, the functions F± entering Eq. ([5]) take the simple form 



F+ 



(tanh(e + eV)/3 ± tanh(e - eV)/3 )/2 



(6) 



where 1/2/3q = To is the reservoir temperature. In the limit of the weak coupling between N and S layers the 
density-of-states in the N metals is supposed to be unperturbed by the proximity effect: A/jv(e) = 1. 

The first terms in the R.H.S. of the kinetic equation (J5]) describe the tunneling of quasiparticles from (to) the normal 
electrodes, whereas the last one is an inelastic collision term. It consists of the electron-electron and electron-phonon 
collision terms (5j„ = S e - e + S e - P h) and is of the order of /±(e)/r e _ e , /±(e)/T e -pfc, where r e _ e and r e _ p /j are the 
electron-electron and electron-phonon inelastic scattering times. Eq.Q contains the DOS, Ng(e,A), depending on 
the energy gap A(V). Therefore Eqs.|l][S|) are coupled, and the problem of finding the distribution function / + is not 
easy. 



III. MODERATE INTERFACE RESISTANCE. HOMOGENEOUS CASE 



We consider first a homogeneous stationary state assuming that all the quantities do not depend on the coordinate 
x. However, the problem of solving Eqs.fl][6]) is too complicated even under this assumption and therefore we consider 
only two limiting cases. 

First, we analyse the limit of a very weak electron-electron and electron-phonon interactions: 

A >> ei tT >> T~_ e , t~_ h , when the collision terms in Eq. ([2]) can be neglected. In other words, the S/N interface 
resistance Rsn is not too high: RsN&d << {-Dr e _ e , Dr e - P h}- This approximation has been used in many papers 
(see, for instance, 0, For a symmetric N/S/N system we use Eq. ([Id]) and obtain from Eq.@ the following 

solutions 



f+=F+, f-=0 



(7) 
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Substituting the function /+ into Eq.([T]) and shifting the energy in the integral by eV in the first term and by — eV 
in the second one, we reduce Eq. |T]) to the following form 

1 = X -Re T tanh ^> ie (8) 

2 J-e» \f{e- eV + ij) 2 - A 2 

Remarkably, the form of Eq. (jSJ) is the same as the one for a superconductor in the presence of an exchange field. In 
other words, the problem of the non-equilibrium superconductor in the N/S/N system is to a great extent equivalent 
to the problem of the equilibrium superconductor with the distribution function /+ = tanh(e/?o) in the presence of 
an "exchange" field eV. 

The latter problem has been attracting a lot of attention since the pioneering works by Larkin and Ovchinnikov 
[34j and by Fulde and Ferrell [33[. For the model with the exchange field, these authors have predicted a new state 
called now LOFF state. They considered a clean superconductor with the exchange field h (or a strong magnetic field) 
acting on spins of electrons. They demonstrated that at zero temperature the energy gap A remained unchanged 
unless the exchange energy h exceeded the value Ao, where Ao is the energy gap at zero temperature in the absence 
of h. 

However, in addition to this solution of the self-consistency equation, there is another, unstable, solution for 
A : A(h) — A \//i 2 //i 2 — 1 with h 2 = Aq/2 and h c < h < Ao. Therefore, in the interval h c < h < Aq there are three 
possible solutions for A : 0, A(h) and Ao (the trivial solution A = always exists). The situation resembles the 
behavior of a non-ideal gas described by the van der Waals equation of state, and thus, one can expect a stratification 
of the electron system. Indeed, as it was shown in Refs. [H, [34J, an inhomogeneous state with the energy gap A(r) 
varying in space turns out to be more favorable than the homogeneous one. 

Using the equivalence between the N/S/N system at a finite voltage V and the equilibrium superconductor in the 
presence of the exchange field one may expect an inhomogeneous LOFF-like state in the system considered here. 
A multi-valued (in a certain region of parameters) dependence of the energy gap A on the applied voltage V and 
damping 7 has recently been found for an N/S/N system in Ref. [23]. The function A(V) had a form similar to that 
determined by Eq. ([8]) . It was established that at some values of the parameters the CVC had the so-called S-shaped 
form. This type of the CVC is well known in bulk superconductors and has already been obtained both theoretically 
(see the review [IS EH) an d experimentally [H,[13] several decades ago. 

Although, one could expect an inhomogeneous state in the N/S/N systems, the authors of Ref. [13] came to the 
conclusion that the homogeneous state had to remain stable. This statement is correct for the superconductor island 
of a small size considered in Ref. [13] but cannot remain valid for large superconducting films sandwiched between 
normal metals. Note also that the authors of Ref. [HJ studied the stability of the system in a short time interval: 
t < h/A (a similar dynamic behavior of the order parameter in a collisionless superconducting system was studied 
earlier in Ref. 43]). On the other hand the instability develops on much longer characteristic times t ~ h/eo. 

In this case, the states corresponding to the part of the CVC with negative differential resistance are unstable 
and the stratification of the current density occurs in the system (see the review [3ljj ) . This structure resembles the 
LOFF coordinate dependence of the order parameter and of other physical quantities. However, in contrast to the 
LOFF state in the equilibrium superconductors with the exchange field, the coordinate dependence cannot be found 
minimizing the free energy because we consider a system out of the equilibrium. 



IV. HIGH INTERFACE RESISTANCE. NON-HOMOGENEOUS CASE 

Although the similarity between the problem involved and the problem of LOFF is most clearly seen in the limit, 
%r >> T~_ e , T~J; ph , one has to solve in this limit complicated integro-differential equations (|1I2[) and this can be done 
only numerically. The problem can be solved analytically in another limit: 

T e-ph « « T-_\ (9) 

In this limit the electron-phonon interaction is very weak and the electron-phonon collision term in Eq. © can be 
neglected. The largest term in Eq. ^ is the electron-electron collision term S e - e {f±}. 

Note that the electron-electron r~_ e and electron-phonon T~_} h scattering rates depend on energy e [35| so that 
Eq.© can be satisfied not at all energies. If the mean free path is not too short (s/T < I, but I < v/3A because 
we consider the dirty limit), these scattering rates are: r~_ e ~ j e - e (T e )c 2 /T 2 and t~ h ~ 7 e _ p / l (T e )e 3 /Tg with 
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FIG. 2: Normalized effective temperature (a) and energy gap (b) as functions of the applied voltage V for 7/A0 = 0.01 and 
different temperatures To = l/(2/3o) (shown in the inset) 

7 e _ e = 7rA e _ e Tg /(8hEp) and j e -ph — ^^e-phT^ / (2hs 2 p F ), where A e _ e , X e - P h are electron-electron (electron-phonon) 
interaction constants, Ep is the Fermi energy, s is the sound velocity (see, for example, [36]). The characteristic energy 
t c h in our case is of the order e c h ~ T e ~ Icj thus, for the Al S-films the inelastic scattering rates are: T~_ e « 10 8 s _1 
and rri^ w 10 6 s -1 [ilij . At temperatures lower than T c , the interval between r~_} e and rjl^, h becomes larger because 
T e-ph nas a stronger dependence on e than r~_} e . However, one has to take into account the dependence of T~_ e and 
t~_ h on A(T). Therefore, if the frequency eo/h = D/(RsNcd), which is determined by the S/N interface resistance 
Rsn, is chosen in the interval 10 6 s _1 — 10 s , the approximation of an effective temperature can be applied to the 
case of these films. 

It is worth noting that the violation of the condition ([9]) leads only to breakdown of the effective temperature 
approximation, but not to disappearance of the effects under consideration. For these effects it is important only to 
have the S-shape CVC, which is realized, for example, also in the case of the frequency eo/h large in comparison with 
r~} e and T~_} ph (see the preceding section and [HI). However the approximation of the effective temperature allows 
one to solve the problem in an analytical form. In the limit of high im pur ity concentration (I < s/T), the inelastic 
electron-phonon scattering rate h depends on energy even stronger |29l. l37l. l38l. [3^| . that is, the window between 

r~} e and t~ h becomes larger and the situation is even more favorable for applicability of our approach. If the 

frequency eo/h is less than T~_ pfl , then the heat absorbed by quasiparticles in the superconductor is released mainly 
to the lattice, but not to the normal electrodes, and our approach is not applicable. However, one can show that in 
this case the CVC of the system also has the S-shape form and all the effects considered here remain qualitatively 
unchanged, although the analytical approach is not possible in this case. 

In the limit determined by Eq. ([9]) the electron-phonon interaction is very weak and the electron-phonon collision 
term in Eq. ^ can be neglected. The largest term in Eq. @ is the electron-electron collision term S e - e {f±}. Due 
to a high rate of the electron-electron collisions one can assume that the distribution function f + has an equilibrium 
form, /+(e) = tanh(e/3), with an effective temperature T = (2/3) -1 depending on the coordinate x. Perhaps for 
the first time, the approximation of the effective temperature has been introduced in the study of noncquilibrium 
superconductors by Parker [4(| (see also [5^, [3l|, SJ). In order to find this temperature, we multiply Eq. ^ by e, 
drop the first term in the L.H.S. and integrate the equation over all energies. This procedure leads to an equation 
describing the energy conservation. The contribution of the collision term S e - e {f+} equals zero because this term 
conserves the total energy, and we obtain the equation for T 
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a ^ f ^->=£<[ W ,f)- S (T)] (10, 



dx 

a—l.r 



where ftf r = (D/ Ri^ad), T = T/Tq is the normalized temperature, and 

M(f) = (f /D) / dyy 2 V + (y, A/3) cosh" 2 ^). 



where the function T> + (y, A/3) is defined in Eq.([4]). 

The functions J(V a , T) and S(T) are the heat source (Joule heat) and drain due to the tunnelling of electrons into 
the N electrodes. They can be written as 

/•OO 

J(V a ,f) = / dyyAf s (y)[tanh(y) - (tanh(y + ) + tanh(y_))/2], (11) 



S(f) = / dyyMs(y)[ttmh(y) - tanh(y/f )], (12) 
Jo 

where y = e/3 , y± = y± eV/3 0l and f = T/T . 

In order to solve Eq. (|10[) we consider for simplicity the symmetric case when V r = — Vi = V. In the homogeneous 
case the dependence of the effective temperature on the applied voltage can be found from the balance equation: 
J{V, T) — S{T). Explicit expressions can be derived analytically in the limits of low and high temperatures. 

In the limit of high temperatures, V f3$ « 1, we come to the following formula 

f-l = (T-T )/T a ^(Vf3 ) 2 (13) 

whereas at low temperatures, V(3o >> 1, \eV — A| >> Tq we obtain 

T-lneV/(A-eV) (14) 

Note that at low temperatures To both the terms J and S are, as they should be, exponentially small but the effective 
temperature T is not. Making these estimations, we assume that 7^0. Thus, it is not obvious that the effective 
temperature T and energy gap A will be multi-valued functions of V as it takes place in another limit considered 
in the previous section. Numerical calculations of the integrals (|11I12[) show that the situation qualitatively remains 
unchanged, i.e. the quantities T and A have three values in a certain interval of voltages V. The dependencies T(V) 
and A(V) are shown in Fig. 2a and 2b. One can see that in a narrow interval of the voltage, the effective temperature 
can have three values (Ti , T2 , T3 ) , and therefore the energy gap A is also a multi- valued function of V", which is 
analogous to the behavior in superconductors with an exchange field h. 
The current through the system is given by the formula 



I= ~ir r ' deJV s (e,A(T))F_(e,V) 
eRb Jo 



(15) 



where F_ = (tanh(e + eV)/3o— tanh(e — eV)Po)/2. The CVC obtained from Eq. (Tl5|) is shown in Fig. 3 for different 
temperatures To (Fig. 3a) and damping 7 (Fig. 3b). 

It is clearly seen from Fig. 3 that the CVC has an "S-shaped" form. This type of the CVC may also occur in 
semiconductors (3l| as a result of the so-called overheating mechanism and in bulk superconductors in the presence of 
a dissipative current [13, EH [H, H^] . This phenomenon is possible for a special dependence of the heat source J(V, T) 
and heat absorption term S(T) on the effective temperature T an d applied voltage V. The "N-shaped" CVC may also 
arise as a result of the overheating mechanism [l2|, HJ, HF 



One can show in the same way as it was done in Ref. 



24,|2J 

3l| that the states corresponding to the part of the CVC 



with a negative differential resistance are unstable and, in the case of a fixed total current [46|, the system is stratified 
into layers with different effective temperatures and current densities. The form of the current filaments can be found 
from Eq. (H 



We introduce a new effective "temperature" d = dT\M(Ti)Ti and a function of this "temperature" W({)) — 
lo \J(Y-> ^(^1) ~ S{T{di))]d , di . Integrating Eq. (fTT)]) over the temperature f and then over d we arrive at the equation 
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FIG. 3: a) The CVC of the system in the homogeneous state for 7/A0 = 0.01 and different temperatures. The part of the 
CVC corresponding to an inhomogeneous state (a broad high temperature domain) is shown by the dashed vertical line, b) 
The CVC in the homogeneous state for different damping 7. The temperature To is choosen so that Aq/3o = 15. 



where l = K r = Ki. 

Eq. (fT6)) is already quite simple and its solution -d describes the dependence of the temperature T on the coordinate 
x along the interface. One can rather easily see that this dependence can be non-trivial. This originates from the fact 
that in a narrow interval of V, where the CVC is a multi- valued function, the function W($, V) has three extrema at 
&k(Tk), where fc = 1,2,3 (see Fig.4a) As a consequence, the solutions of Eq. HU), l Q {dd/dx) = V^y/W - W{^,V), 
have different forms: solitons (instantons), oscillatory temperature distributions or domain walls. 

In Fig. 4b we show qualitatively phase trajectories illustrating this conclusion. At a certain voltage Vo the function 
W{d) has the same values at the maxima: W($i) = W{$s). This means that the voltage Vo satisfies the condition: 
Jtf? [J{Vo,T(0) - S(T(9))]de = 0. The trajectory connecting these maxima (separatrix) corresponds to a solution like a 

domain wall: l^dd/dx) = \f2y/W{d\, Vq) — W("d, Vo). The spatial distribution of the effective temperature T related 
to this trajectory is a decay of T from the value T3 to T\ over a length of the order of 1$. The trajectories close to 
the separatrix describe the domain structure (see Fig. 4c). This structure is analogous to the one in a superconductor 
with an exchange field in the LOFF state. 

At the same time, the inhomogeneous state in the SN systems out of the equilibrium differs from the LOFF state. 
In the latter case the solution of the self-consistency equation should correspond to the minimum of the free energy 
whereas in our case the choice of the solution is dictated by the bias current lb and by the stability against small 
perturbations. One can show that the most stable solution has the minimal number of zeros of the function dT(x)/dx 

This fact suggests the following scenario: 

a) With increasing the bias current lb the dependence /(V) follows the CVC for the homogeneous case. 

b) At lb > I m the voltage V jumps to a lower value close to Vq and remains almost constant with increasing lb (see 
Fig. 3a and [3l|). 

c) When the current h exceeds the value I3, the voltage V increases along the CVC for the homogeneous case. If 
lb decreases from ^3, the voltage decreases and a jump to a higher voltage close to Vo occurs at lb < Im (a hysteretic 
behavior). 

It would be interesting to study experimentally the inhomogeneous distribution of the effective temperature T(x), 
current and the gap A(x). This could be done, for example, by applying the point-contact spectroscopy to N/S systems 
(see a possible setup in Fig. lb). The distribution function / + in these systems has the same form as in N/S/N systems. 



(l/2)ll(d$/dx) 2 =W Q - W(0, V) 



(16) 
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FIG. 4: a) The potential W(6) for different V: V = V (solid line); V < Vo (dashed line) and V > V (dash-dotted line), b) 
Phase trajectories (T, dT/dx). Solid line shows a separatrix (a domain wall in T{x) distribution); dashed and dotted lines show 
domain-like structure of T(x) distribution; c) Spatial dependence of the effective temperature corresponding to a domain-like 
structure (arbitrary units) 



Therefore, by measuring the spatial distribution A(x), one can get information about the inhomogeneous states in 
the N/S systems. 

It is interesting to note that the ideas that the symmetry breaking of the LOFF type may occur not only in the 
superconductors but also in the quantum chromo-dynamics, astrophysics and in cold gases (see the review [48| and 
references therein) have been put forward recently. Many theoretical papers have been published since the pioneering 
works [lil, [34[ , in which different aspects of the LOFF state in superconductors and conditions necessary for this state 
were analyzed (see Refs. [H, [4{| and references therein). As to the experimental observation of the inhomogeneous 
states predicted, the situation is not as clear. Although some observations in low dimensional superconductors and 
superconductors with heavy fermions can be interpreted in terms of the LOFF states, there are no convincing evidences 
in favor of such states in ordinary s-wave superconductors. 

So, the investigation of non-equilibrium superconductors, for example superconductors in tunnel N/S/N systems, 
can provide one more possibility to observe the LOFF state experimentally. Such experiments may be useful for 
understanding the LOFF states and the conditions under which they can be realized. We mention here another inter- 
esting example of the analogy between a non-equilibrium superconducting system and equilibrium superconducting 
system with an exchange field, namely, the 7r— state may arise not only in S/F/S junctions [5(| HH H3] , but also in 
S/N/S junctions with a non-equilibrium distribution function (l5l. Hfl UtI Il8j (here F denotes a ferromagnetic layer). 



V. CONCLUSIONS 



In conclusion, using the approximation of an effective temperature of quasi-particles, we have studied non- 
equilibrium states in a tunnel NSN structure with a low barrier transparency. It is found that for certain values 
of parameters, the CVC of the system may have an S-shaped form. The uniform state corresponding to the part of 
the CVC with negative differential conductance is unstable, and therefore a nonuniform current I{x) and temperature 
T(x) distribution is established in the system with a fixed total current. We discuss the analogy with the nonuniform 
LOFF states in equilibrium superconductors and the possibilities of experimental observation of the nonuniform states 
in SN structures. 
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